AD-A230  560 


Conversion  of  Electrostatic  and  Electromagnetic  Waves  in  a 
Plasma  at  the  Peak  of  a  Parabolic  Density  Profile 


D.E.  Hinkel-Lipsker,  B.D.  Fried,  and  G.J.  Morales 

PPG-1338  December  1990 


JAN  03  1991 


University  of  California  at  Los  Angeles 

Department  of  Physics 

Los  Angeles,  CA  90024-1547 


\ 


t  ' 


1 


Conversion  of  Electrostatic  and 
Electromagnetic  Waves  in  a  Plasma  at  the  Peak 
of  a  Parabolic  Density  Profile 

D.E.  Hinkel-Lipsker,  B.D.  Fried,  and  GJ.  Morales 
Department  of  Physics,  University  of  California,  Los  Angeles 
Los  Angeles,  California  90024-1547 

Abstract 

Analytic  expressions  for  the  reflection,  transmission,  and  mode  conversion  co¬ 
efficients  for  electromagnetic  and  Langmuir  waves  in  an  unmagnetized  plasma 
with  a  parabolic  density  profile  are  found  for  both  the  “direct”  problem  (incident 
electromagnetic  wave)  and  the  “inverse”  problem  (incident  Langmuir  wave).  In 
contrast  to  the  linear  profile  problem,  the  absorption  depends  explicitly  on  the 
value  of  the  collision  frequency  (cold  model)  or  temperature  (warm  plasma),  but 
a  transformation  of  parameters  relates  the  results  for  these  two  limits. 


PACS  Numbers:  52.40Db,  52.35Lv,  94.20Bb 


The  conversion  of  electromagnetic  waves  into  electrostatic  modes  in  an  inho¬ 
mogeneous,  unmagnetized  plasma  is  a  topic  of  broad  interest  which  has  impor¬ 
tant  consequences  for  laboratory  and  space-plasma  studies.  Most  of  the  analysis 
of  this  problem  has  concentrated  on  the  simpler  case  of  a  plasma  with  a  linear 
density  profile1-7;  in  fact,  only  recently6  has  an  analytic  solution  been  obtained 
for  this  case  by  a  scale-separation  technique.  The  present  study  uses  this  tech¬ 
nique  to  solve  the  more  technically  demanding  parabolic  density  profile  problem, 
a  physically  important  situation  of  current  interest  to  ionospheric  heating  at  the 
peak  of  the  F-layer*,  Langmuir  turbulence  studies  in  basic  laboratory  devices9, 
beat  wave  acceleration  schemes,  and  some  aspects  of  laser  plasma  interactions. 

Previous  analytical  studies  of  this  problem  have  considered  the  structure 
of  the  electromagnetic  wave,  as  done  by  Baftos10  using  Langer’s  method,  and 
by  Lundborg  and  Thidi 11,12  through  an  analytic  uniform-approximation  method. 
However,  these  investigations  have  not  yielded  the  structure  of  the  electrostatic 
fields  excited  at  plasma  resonance,  and  have  not  predicted  the  experimentally  im¬ 
portant  mode  conversion  coefficients.  The  present  study  yields  relatively  simple 
analytic  expressions  for  these  previously  unknown  quantities  and  simultaneously 
solves  die  inverse  mode  conversion  problem,  a  situation  of  diagnostic  interest 
since  it  permits  the  sampling  of  microscopic  plasma  phenomena  at  remote  loca¬ 
tions. 

In  the  direct  problem,  an  electromagnetic  wave  of  frequency  uj  is  obliquely 
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incident  on  the  plasma  at  an  angle  0  relative  to  the  density  gradient,  with  u>  cho¬ 
sen  equal  to  the  electron  plasma  frequency  wp  at  the  peak  of  the  parabolic  density 
profile  (Fig.  1).  This  wave  encounters  a  cut-off  at  the  layer  wp(z)  =  u>cos0,  where 
some  of  the  wave  is  reflected,  while  the  remainder  tunnels  to  the  plasma  reso¬ 
nance  where  u >  =  o>p(0).  In  this  region,  a  portion  of  the  incident  wave  energy  is 
transformed  into  plasma  oscillations  for  the  cold,  collisional  problem,  and,  for  a 
warm  plasma,  into  Langmuir  waves  having  cut-offs  at  u>p(z)  =  u>(l  -  /32sin20)'/2, 
where  (32  =  3Te/me(P’.  The  remaining  energy  tunnels  out  to  the  electromagnetic 
cut-off,  where  it  becomes  the  transmitted  fraction  of  the  electromagnetic  wave. 
Fig.  1  depicts  this  process  for  the  warm,  direct  problem. 

We  also  consider  the  warm,  inverse  problem,  where  an  incident  (leftgoing  in 
Fig.  1)  electrostatic  wave  in  the  plasma  gives  rise  to  two  outgoing  electrostatic 
waves  [a  rightgoing  wave  resulting  from  the  reflection  at  the  electrostatic  cut-off 
and  a  leftgoing  wave  transmitted  through  the  evanescent  layer  around  z  -  0]  and 
two  obliquely  propagating,  outgoing  electromagnetic  waves  produced  by  mode 
conversion  in  the  region  around  z  =  0. 

In  general,  the  electrons  are  treated  as  a  warm,  collisional  fluid  with  an 
unperturbed  density  n0(z)  =  n(l  -  z2/L2)  and  collision  frequency  v.  The  ions 
are  static,  and  quasineutrality  is  assumed.  The  fluid  equations  are  linearized  by 

setting  n(F,  t)  -  n„  +  {ni(z)exp[tkxz  -u;f]  +  c.c.},  etc.  We  assume  that  there  are 

— ♦  — # 

no  zero  order  fields  or  flows  ( Ea  =  Ba  =  u0  =  0).  The  fluid  equations,  together 
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with  Maxwell’s  equations  then  yield  a  set  of  equations6,7  which  couple  Et  and 
By 

This  set  of  equations  can  be  expressed  in  terms  of  scaled  coordinates  and 
wavenumbers  appropriate  for  the  parabolic  profile  problem:  Z  and  N  for  the 
electromagnetic  mode,  and  (  and  K  for  the  electrostatic  mode,  where 

C  =  k'J2z/(l3L)'12  =  Z/f3{'2,  (1) 

and 

K  =  f3'/2(kxL)/(k0L)l/2  =  (Sl/2N  ,  (2) 

with  k0  =  u>/c . 

When  the  frequency  matches  the  plasma  frequency  at  the  peak  of  the  profile, 
the  scaled  set  of  equations  have  the  form 

K  +  «2  +  ivo/0  -  K2)Et  =  bB(Z),  (3) 

p"  +  «2  +  iv0/(3  -  K2)p  =  -2 CEJa,  (4) 

d2B/dZ2  -  2 Z(Z2  +  iv„  -  qr'dB/dZ  +  (Z2  +  iu„  -  q)B  =  5,  (5) 

where  koP  =  V  E,  B  =  iBv ,  b  =  iNa(l  -  (32)0~V2,  a  =  (k JJ7?'2, 
=  jb0Xi//w,  5  =  S{Z)= -HNZQ{Z2+iv0-q),  with  Q  =  (1  ~f32)P-a~' E'z, 
q  =  N2  =  (k^Lj2 /(k0L),  and  the  prime  denotes  differentiation  with  respect  to 
Eqs.  (3)  -  (5)  describe  wave  propagation  and  conversion  in  a  plasma  with 
a  parabolic  profile.  The  electrostatic  modes  vary  rapidly,  i.e.,  on  the  short 
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scale  (0.  whereas  the  electromagnetic  modes  vary  much  more  slowly,  i.e.,  on 
the  longer  scale  Z.  The  coupling  of  the  two  modes  depends  on  the  angle 
of  incidence  through  N.  The  larger  the  angle  of  incidence,  the  stronger  the 
coupling;  however,  as  N  increases,  so  does  the  tunnelling  distance  to  the  mode 
conversion  point,  and  the  competition  between  these  two  phenomena  causes  a 
peak  in  the  mode  conversion  coefficient  for  a  finite  angle  of  incidence. 

Cold,  Collisional  Problem.  For  this  case,  we  let  f32  =  0  in  the  set  of 
equations  (3)  -  (5),  but  keep  v0  J  0.  In  this  limit,  equations  (3)  and  (4)  yield 
algebraic  relations  for  p  and  Et  in  terms  of  B.  When  these  expressions  are 
substituted  into  the  righthand  side  of  (3),  we  obtain  a  differential  equation  for 
B  whose  homogeneous  solutions  are  derivatives  of  parabolic  cylinder  functions. 
This  equation  is  solved7  using  a  Green’s  function  technique,  together  with  the 
approximation  that  on  the  righthand  side  of  (3),  which  can  be  shown  to  be 
proportional  to  B(Z),  we  replace  5(Z)  by  5(0).  We  then  self-consistendy 
determine  5(0)  as  well  as  5(Z)  for  Z  »  l  and  Z  <  —  1,  from  which  the 
reflection  and  transmission  coefficients  can  be  found.  We  also  compute  the 
mode  conversion  coefficient  |?7|2  (the  ratio  of  power  collisionally  dissipated  in 
plasma  oscillations  to  the  incident  power  in  the  electromagnetic  wave)  which 
satisfies  the  energy  conservation  equation 

1111  =  1  -  |Ji|J  -  |7f.  (6) 
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The  final  results  can  be  most  succinctly  expressed  in  the  form 

i4jrr 


\T\2  =  (l  +  e?’r' 
\v\2  = 


i/^irp  +  .Yi2’ 
2u\t2\T\2  Im[(i  +  e^/2)X] 


(7) 


(8) 


\v0/2\T\2  +  X\*  l+e^ 
where  X  =  2g7r2exp[97r/4  +  3?ri/4]/[t  +  exp(<77r/2)],  |T|  =  |F[(1  +i^)/4]|  is  the 

gamma  function13,  and  |JR|2  is  determined  from  (6).  The  converted  energy  from 
just  one  side  of  the  mode  conversion  layer  is  half  of  the  amount  given  by  (8). 
Note  that  \q\2  vanishes  as  va  — *  0,  in  contrast  to  the  linear  profile  problem  where 
I77I2  has  a  non-zero  value  in  the  cold,  collisionless  limit.  In  the  limit  q  — ♦  0, 
|t;|2  *  0,  and  |f?|2  =  |T|2  =  0.5,  in  agreement  with  the  earlier  work  of  Banos.10 

Fig.  2  depicts  the  dependence  of  the  energy  flux  coefficients  on  the  parameter 
q  =  kaL  sin20  for  v0  -  .01.  For  this  case,  the  peak  of  |r;|2  occurs  at  q  = 
qa  =2  0.05,  which  is  smaller  than  for  the  linear  profile  problem67,  where  q\  = 
(koL)2/3sin20  0.5.  [It  follows  from  (8)  that  for  small  q,  \q\2  is  proportional 

to  q/\q  +  a,y2|2,  where  a  is  of  order  10  and  is  independent  of  q,  so  jr; |2  peaks 
at  a  value  of  q  as  For  larger  v>a,  the  peak  of  \q\2  moves  to  higher  q,  as 

does  the  minimum  of  |iZ|2,  and  the  tail  on  the  mode  conversion  coefficient  falls 
off  more  slowly.  As  q  increases,  the  tunnelling  distance  for  the  electromagnetic 
wave  becomes  large.  This  thick  evanescent  layer,  which  causes  most  of  the 
incident  energy  to  be  reflected,  accounts  for  the  fall-off  of  |T|2  and  |?7|2  and  die 
large  value  of  |f?|2  (approaching  1)  seen  in  Fig.  2  for  q  >  1. 
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Warm,  Direct  Problem.  In  this  case  we  set  u0  =  0  and  keep  (3  ^  0.  In 
(3),  the  field  5,  is  on  the  short  electrostatic  scale  (()  whereas  its  source  on  the 
righthand  side,  which  is  proportional  to  5(Z),  is  on  the  longer  electromagnetic 
scale.  We  utilize  this  disparity  in  scales  to  replace  5(Z)  by  5(0)  on  the  righthand 
side  of  (3).  Equation  (3)  is  then  just  a  parabolic  cylinder  equation  driven  by  a 
constant  source.  The  solutions  are  determined  by  using  the  Green’s  function  for 
outgoing  waves. 

Given  this  solution  to  Ez  [in  terms  of  5(0)]  and  the  boundary  conditions 
for  the  direct  problem,  it  immediately  follows  that  p(0  =  a~'5'.  The  variables 
Et  and  p  thus  have  functional  forms  expressible  in  terms  of  parabolic  cylinder 
functions  with  argument  (  and  their  derivatives,  with  amplitudes  proportional  to 
the  quantity  5(0). 

In  order  to  determine  5(0),  we  substitute  p{Q  and  5,(0  into  the  source 
term  on  the  righthand  side  of  (5).  This  equation  has  as  homogeneous  solutions 
the  derivatives  of  parabolic  cylinder  functions,  and  these  are  used  to  form  the 
Green’s  function  for  B(Z)  satisfying  the  boundary  conditions  that  5(Z)  con¬ 
sist  of  leftgoing  (incident)  and  rightgoing  (reflected)  waves  for  Z  >  1,  and  a 
leftgoing  (transmitted)  wave  for  Z  <  -1.  We  then  make  the  approximation 
that  the  parabolic  cylinder  functions  with  argument  Z,  which  are  multiplied  by 
5,(0  and  p(Q  in  the  Green’s  function  expression  for  5,  can  be  replaced  by 
their  Z  =  0  values  during  integration  on  the  rapidly  varying  (  scale. 
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From  the  resulting  expression  for  B(Z),  we  seif-  consistently  determine  5(0) 
[by  evaluating  B(Z  =  0)]  and  then  calculate  B(Z  >  1)  and  5(Z  <  -1),  which 
give  the  reflection  and  transmission  amplitudes,  respectively.  The  results  can 
most  easily  be  expressed  as  follows.  Each  of  the  quantities  |T|2,  \R\2,  and 
|r/|2  for  the  warm,  direct  problem  is  equal  to  the  corresponding  quantity  for  the 
cold,  collisional  problem  with  the  substitutions  v0  — -  0  and  X  — ►  inM 2,  where 
M2  =  ^7eiexp[(l  +  /3)qir/4][fl™  7(0<£<]/{[t  +  exp(gtr/2)][i  +  cxpiPqir /2)]}, 
and  7«  =  m  +  i<7/?)/4]/r[(l  -  i<7/3)/4].  The  function  I(Q,  which  represents 
the  solution  to  the  driven  parabolic  cylinder  equation  with  a  constant  source,  is 
analogous  to  the  function7  Gi  +  iAi,  which  is  an  outgoing  solution  to  the  driven 
Airy  equation,  encountered  in  the  linear  density  profile  problem. 

Fig.  3  exhibits  the  warm  energy  flux  coefficients  as  a  function  of  q  for 
electron  temperature  f32  =  10-4.  Since  M2  and  X  have  essentially  the  same  q 
dependence  as  long  as  f3q  <  1,  it  is  not  surprising  that  these  curves  are  similar 
to  those  of  Fig.  2  for  the  cold,  collisional  case.  In  both  the  cold  and  warm 
cases,  the  small  q  form  of  jr/l2  is  abq/\a  +  ibq |2,  where  a  and  b  are  independent 
of  q,  so  the  maximum  value  of  \tj\2  is  0.5. 

Warm,  Inverse  Problem,  Solution  of  the  inverse  problem  proceeds  in  a 
manner  analogous  to  that  of  the  direct  problem,  with  one  technical  complication 
which  is  a  consequence  of  the  difference  in  boundary  conditions  for  the  direct  and 
inverse  problems.  In  the  solution  to  the  p  equation,  (4),  the  boundary  conditions 
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for  the  inverse  problem  permit  a  homogeneous  solution  as  well  as  the  particular 
solution.  The  coefficient  p„  of  this  homogeneous  solution  constitutes  a  second 
unknown  constant,  in  addition  to  B(0).  We  determine  these  by  evaluating  the 
Green’s  function  expression  for  B(Z)  at  Z  =  0,as  in  the  direct  problem,  and 
also  using  the  2  component  of  Ampere’s  Law  evaluated  at  Z  =  0.  We  find  that 
the  effect  of  p„  is  simply  a  correction  of  order  0{/2  to  the  constant  2?(0),  so  we 
can  neglect  the  p0  term. 

The  energy  flux  coefficients  for  the  inverse  problem,  shown  in  Fig.  3  as 
a  function  of  q  for  02  =  104,  are  determined  in  a  manner  identical  to  that 
for  the  direct  problem.  Only  a  single  mode  conversion  coefficient  is  plotted, 
since  we  have  shown  by  explicit  calculation  that  \rj§  =  |r/|£  (where  the  I  and 
D  subscripts  refer  to  the  inverse  and  direct  cases,  respectively),  a  result  which 
follows  from  time  reversal  symmetry.  Fig.  3  also  reveals  an  additional  near 
symmetry:  |/?|J  %  |T|d  and  \T\j  ~  |i2|J.  We  have  analytically  demonstrated 
that  these  equalities  are  valid  to  within  terms  of  order  0l/2,  but  have  not  been 
able  to  show  exact  equality. 

From  the  present  results  it  is  clear  that  the  degeneracy  between  cold  and 
warm  plasma  behavior  found  in  the  linear  profile  problem  is  not  present  for  the 
parabolic  profile.  In  the  latter  case,  all  mode  conversion  coefficients  depend 
explicitly  on  v„  and  0.  However,  it  is  possible  to  relate  the  warm  absorption 
coefficients  to  the  collisional  results  through  a  transformation  of  parameters.  For 
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non-linear  studies,  of  relevance  to  ionospheric  and  laboratory  experiments,  it  is 
worth  emphasizing  that  for  the  same  input  power,  the  experimentally  measurable 
ratio  of  peak  electric  fields  is 

1-£z1pAR  _  _ 18.6 _ ^PAR  Q 

I^Iun  (/?ff)'/2(13.i/?+.05)‘/2ZuN  ’  } 

where  Lus  is  the  scale  length  of  the  linear  density  profile  and  Lpar  is  the 
curvature  of  the  parabolic  density  profile. 

The  present  study  clearly  illustrates  the  usefulness  of  the  scale-separation 
technique.  It  allows  an  analytic  solution  of  the  direct  and  inverse  conversion 
problems  for  both  linear  and  parabolic  density  profiles  from  a  unified  point  of 
view.  We  believe  that  the  method  may  be  of  value  in  similar  situations  involving 
magnetized  plasmas. 

This  work  was  supported  by  the  Office  of  Naval  Research,  U.S.DOE,  and 
the  UCLA  Dissertation  Year  Fellowship. 
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Figure  Captions 


Figure  1.  Schematic  of  the  direct  problem.  An  incident  electromagnetic  wave 
(EMj)  gives  rise  to  a  reflected  wave  (EMr),  a  transmitted  wave  (EM-r),  and 
outgoing  Langmuir  waves  (ESout)-  The  inner  and  outer  dashed  lines  represent 
the  electrostatic  and  electromagnetic  cut-offs,  respectively. 

Figure  2.  Cold,  collisional  energy  flux  coefficients  as  a  function  of  q  =  k0Ls'm29 
for  ua  =  .01. 

Figure  3.  Warm,  energy  flux  coefficients  as  a  function  of  q  ~  JfcoZsin20  for  (32  = 

io-4. 
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